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Three elastic phases of covalent networks, (I) floppy, (II) isostatically rigid 
and (III) stressed-rigid have now been identified in glasses at specific degrees 
of cross-linking (or chemical composition) both in theory and experiments. 
Here we use size-increasing cluster combinatorics and constraint counting al- 
gorithms to study analytically possible consequences of self-organization. In 
the presence of small rings that can be locally I, II or III, we obtain two tran- 
sitions instead of the previously reported single percolative transition at the 
mean coordination number f = 2.4, one from a floppy to an isostatic rigid 
phase, and a second one from an isostatic to a stressed rigid phase. The width 
of the intermediate phase Af and the order of the phase transitions depend 
on the nature of medium range order (relative ring fractions). We compare 
the results to the Group IV chalcogenides, such as Ge-Se and Si-Se, for which 
evidence of an intermediate phase has been obtained, and for which estimates 
of ring fractions can be made from structures of high T crystalline phases. 

Pacs: 61.20N-81.20P 
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I. INTRODUCTION 



A covalently bonded amorphous network progressively stiffens as its connectivity or mean 
coordination number f increases^'^. The increase of connectivity can be achieved by adding 
cross-hnking elements (such as As or Ge) to a starting chain network of S or Se. From a 
mechanical viewpoint, two- fold coordinated single bond chain networks are floppy because 
the number of nearest neighbor (r=2) bonding constraints per atom is less than three, the 
degrees of freedom^. On the other hand, in the case of networks consisting only of tetrahedral 
units (such as amorphous four-fold coordinated silicon), the network is intrinsically rigid. 

In single bond random networks, these simple observations have been described with 
success with a mean- field theory based on Maxwell constraint counting^. For r-fold co- 
ordinateded atoms, the enumeration of atomic constraints — r/2 and — 2r — 3, 
respectively, due to bond stretching and bond-bending forces has shown that the number 
of zero-frequency (floppy) modes per atoms actually vanishes when the mean coordination 
number^'^ of the network increases to the magic number of 2.4. At this point, the network 
sits at a mechanically critical point where the number of constraints per atom ric = n" + 
equals the number of degrees of freedom per atom. On the one hand, the exhaustion of 
all the degrees of freedom means that the network efficiently fills space. At the same time, 
because there are no excess constraints, the network can be thought of as globally isostatic. 

Mean-field theory predicts onset of rigidity for f > 2.4. Numerous experiments^'^^'^^ 
have confirmed these simple predictions, especially in glass science where bulk chalcogenide 
glasses have been used as a benchmark to check these elegant ideas. Threshold behavior 
has been detected in structural^, vibrational^, thermal^" and electronic^^ properties when 
f approaches 2.4. Applications of rigidity theory have also been reported in various fields 
such as granular matter, biology and computational science^^. However, recently it has 
been shown from Raman scattering and from phase-dependent measurements of the kinet- 
ics of the glass transition that two transitions^^ at fd and fc2 appear when the network 
stiffens. This suggests that the mean-field constraint counting alone (leading to the single 
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percolative transition) fails to describe completely network changes. These two transitions 
define an intermediate phase in which the connected structure continues to be stress-free^^'^^ 
(isostatically rigid). 

In this paper, we present a simple way to go beyond the mean-field description of rigidity 
and include local stress corrections. This is achieved by performing Maxwell constraint 
counting on size increasing structures, starting from the short-range level corresponding to 
the previous mean-field approach. One main advantage is that medium range order effects 
such as small rings can be taken into account in this construction. It appears from this 
analysis that these small rings mostly determine the nature of the intermediate phase and the 
values of the critical coordination numbers fd and fc2, hence the width of the intermediate 
phase Af — fc2 — fd- To apply this construction, we choose the simplest case that can be 
built up, and which has received considerable attention in the context of rigidity, namely 
single-bonded Group IV chalcogenide glasses of the form B^Ai^x with coordination numbers 

= 2 and = 4 defining the mean coordination number f — 2 + 2x. We have used size- 
increasing cluster approximations (SICA) to construct these size-increasing structures and 
medium range order (MRO) on which we have realized constraint counting. As hoped, the 
analysis reveals two transitions: a first one, at which the number of fioppy modes vanishes, 
is closely related to what has been previously obtained in the mean-field approach, and a 
second one (a "stress transition"), beyond which stress in the entire structure can no longer 
be avoided. The second transition can be obtained only under certain conditions which we 
detail below. We show that the orders of these phase transitions are different and depend 
on the fraction of ring structures. In between the two transitions, one can define an almost 
stress-free network structure for which the fraction of isostatic clusters can be computed. 

The paper is organized as follows. In Section II, we show how to construct from small 
molecules size increasing clusters using SICA and perform Maxwell constraint counting on 
them. Section III is devoted to the results obtained from the construction, and the change 
in structure and energy with increasing connectivity. We discuss the results obtained and 
compare them with chalcogenide glasses in Section IV and finally, we extend the approach 
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to fast ionic conducting glasses in Section V. 



II. SIZE INCREASING CLUSTER APPROXIMATIONS 

A. Construction 

In this section, we describe size increasing cluster approximations (SICA). This approach 
has been first introduced to describe the ring statistics and the intermediate range order in 
amorphous semi-conductors^^ such as B2O3 but other applications have been considered 
such as the high temperature formation of fuUerenes^^ and the cell distribution in quasi- 
crystals^°. In principle, any structural quantity that is computed when one increases the 
size of a given structure (or a starting network) converges to its "true" value if the size 
becomes almost infinite. In practice, one hopes that the convergence is rapid enough to give 
reasonable values for medium-sized clusters, yielding information about MRO structures. 
In quasi-two-dimensional v — B2O3 clusters having ten boron atoms allows one to obtain 
a fraction of boron atoms trapped in boroxol rings which is in very fair agreement with 
experiments^, but larger clusters may be needed for three-dimensional networks. 

The basic level of the SICA construction is the restricted mean-field approximation 
where the probability of the short range order structure is derived from the macroscopic 
concentration, assuming that the cations and anions alternate in the network (chemical 
ordering). This basic level is denoted by / = 1. Then, we construct I = 2, I = 3 etc. and 
compute the corresponding probabilities in the Canonical Ensemble with particular energy 
levels Ei corresponding to bond creation between the {I — 1) short range order molecules 
which are used as building blocks from step (/ = 1) (corresponding in GcxSei-x also to the 
reported mean-field approach^) to arbitrary 1. The construction is supposed to be realized at 
the formation of the network, when T equals the Active temperature Tf so that Boltzmann 
factors of the form = exp[—Ei/Tf] will be involved in the probabilities^^ . Since we 
expect to relate the width of the intermediate phase to the ring fraction, we will restrict 
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our present study to Group IV chalcogenides of the form SixSei-x- For the latter, there 
is strong evidence that at the stoichiometric concentration x=0.33 a substantial amount of 
edge-sharing SiSe4^/2 tetrahedra^^'^"^ can be found. 

In order to study Group IV chalcogenides, we select basic units such as the A2 (i.e. 862) 
chain fragment and the stoichiometric ^^44/2 molecule (e.g. GeSe4/2 which is the majority 
local structure at a; = 0.333). These basic units have respective probabilities 1 —p and 
p — 2x/{l — x), X being the macroscopic concentration of the Group IV atoms. The energy 
levels are defined as follows. We associate the creation of a chain-like A2 — A2 structure with 
an energy gain of Ei, A2 — BA2 bondings with an energy gain of E2 and corner-sharing (CS) 
and edge-sharing (ES) -BA4/2 tetrahedra or any ring structure respectively with £'3 and £'4. 
The energy £4 will be used to change the fraction of edge-sharing relative to corner-sharing 
tetrahedra. The probabilities of the different clusters have statistical weights g{Ei) which 
can be regarded as the degeneracy of the corresponding energy gain and correspond to the 
number of equivalent ways a cluster can be constructed. For the step 1 — 2, four different 
clusters can be obtained (see fig. 1) and their probabilities are given as follows: 

^ 4(1 . y 

4(1 -p)2ei + 16p(l-p)e2+p2(16e3 + 7264) ^ ' ' 

P2 = mi~p)e2 

4(1 -p)2ei + 16p(l-p)e2 + p2 (1663 + 7264) 

P3 = (2.3) 

4(1 -p)26i + 16p(l-p)62+p2(1663 + 7264) ^ ^ 

- 72p'^e4 

~ 4(1 - p)26i + 16p(l - p)62 + p2(l663 + 7264) ^ ' ' 

out of which the concentration of B atoms x^"^^ can be computed: 

^(2) ^ P2 + 2(P3+P4) ^2 5^ 

4-p2-2(p3-P4) ^ ^ 

Due to the initial choice of the basic units, the energy E2 will mostly determine the proba- 
bility of isostatic clusters since the related Boltzmann factor 62 is involved in the probability 



(2.2) of creating the isostatic BA4 cluster (a A2 — BA4/2 bonding). This means that if we 
choose to have E2 <^ Ei, E3, E4, the network will be mainly isostatic in the range of interest. 

For the next steps, care has to be taken in order to count the possible isomers produced 
from different pathways (e.g. in fig. 1, the six-membered ring with two B atoms can be 
produced out of p2 and ps). More generally, increasing steps will lead to clusters with 
stoichiometry GenSe2i and probability proportional to p"(l — p)^~" with n = 0..1. The 
corresponding statistical weights depend much more on the way the clusters are created 
and have therefore no general formula depending on n or I. However, for the pure edge- 
sharing tetrahedra chain, it can be easily checked that its probability is proportional to 
72 X 24'~^p'e4~^. Another simple example is provided by the Se chain whose probability is 

It is obvious that all the cluster probabilites will depend only on two parameters (i.e. the 
factors ei/e2 and 63/62) and eventually 64/62 if one considers the possibility of edge-sharing 
(ES) tetrahedra or rings. One of the two factors can be made composition dependent since 
a conservation law for the concentration of B atoms x^''^ can be written at any step / of the 
SIC A construction^^: 

x^') = X (2.6) 

This means that either the fictive temperature Tf or the energies Ei depends^ ^ on x but 
here only the ei{x) (or 6j(f)) dependence is relevant for our purpose. 

The construction has been realized up to the step / = 4 which already creates clusters of 
MRO size. 

B. Maxwell cluster constraint counting 

On each cluster one can count Maxwell constraints by enumeration of bond-bending and 
bond-stretching constraints and calculation of the corresponding expressions of and . 
Of particular importance are the structures containing a ring (see Fig. 1), because special 
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care has to be taken to avoid the counting of redundant constraints^. To illustrate this, 
let us consider an isolated triangle (i.e. a three-membered ring) having a 2-fold atom at 
each of its vertices. Of course, this triangle can be completely defined by three independent 
variables (e.g. two lengths and one angle). Performing constraint counting on the atoms 
will give three bond-stretching and three bond-bending constraints, yielding three redundant 
constraints. This means that for a three-membered ring, one has to remove three constraints 
from the global counting. For the four-membered ring, this correction is of two constraints, 
and for a five-membered ring, of one. 

For each step /, we have computed the total number of constraints n|,: 

n« = (2.7) 

where nc{i) and iVj are respectively the number of constraints and the number of atoms of 
the cluster of size I with probability pi. Hi is the total number of clusters of size I. At step 
Z = 2, it is easy to check that: 

^(2) ^ 4pi + 15p2 + 22p3 + 20p4 ^2 8) 

2pi + 5p2 + 6(p3+P4) 

We have determined either ei/e2 or 63/62 by solving equ. (2.6), and once these factors 
become composition dependent, it is possible to compute the probabilities pi as a function 
of composition and find for which concentration x (or which mean coordination number f) 
the system reaches optimal glass formation where the number of floppy modes fi — 2> — n^l^ 
vanishes. 

III. RESULTS 

A. Structural properties 

In this section, we consider the solutions of the SICA construction under various struc- 
tural possibilities. 
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The simplest case is random bonding, which is obtained when the cluster probabilities 
Pi are only given by their statistical weights g{Ei). This would for example reduce the 
probability p4 in equ. (2.4) to 

with p = 2(f — 2)/(4 — f). A single solution is obtained for the glass optimum point defined 
by the vanishing of the number of floppy modes / = at all SICA steps, in the mean 
coordination number range [2.231,2.275], shghtly lower than the usual mean-field value of 
2.4. Since there is only one solution, there is no intermediate phase in the case of random 
bonding. 

Self-organization of the network can be obtained by starting from a floppy cluster of size 
/ (e.g. a chain-like structure made of a majority of A atoms), and allowing the agglomeration 
of a new basic unit onto this cluster to generate the cluster of size l+l only if the creation of a 
stressed rigid region can be avoided on this new cluster. This happens when two BA4/2 basic 
units are joined together on a given cluster. With this rather simple rule, upon increasing 
f one accumulates isostatic rigid regions on the size increasing clusters because BA4/2 units 
are only accepted in A2 — BA4/2 isostatic bondings with energy E2. On the opposite side, 
one can start at high concentration, close to the mean coordination number of f = 2.67 and 
follow the same procedure but in opposite way, i.e. with adding A atoms, one allows only 
bondings which lead to isostatic rigid or stressed rigid regions, excluding systematically the 
possibility of floppy A2 — A2 bondings. 

In the case of self-organized clusters, the simplest case to be studied is the case of dendritic 
clusters, where no rings are allowed (achieved by setting 64 to zero). For an inflnite size I, 
this would recover the results from Bethe lattice solutions or Random Bond Models^^ for 
which rings are also excluded in the thermodynamic limit^^. A single transition for even I 
steps at exactly the mean-fleld value f = 2.4 is obtained whereas for the step I — 3, there is 
a sharp intermediate phase deflned by / = (still at f = 2.4) and the vanishing of floppy 
regions (i.e. ei/e2 is zero) at f = 2.382(6). Once the probabilities of floppy, isostatic rigid 



and stressed rigid clusters as a function of the mean coordination number are computed, it 
appears that the network is entirely isostatic at the point where f — (solid line, fig. 2). 
There the number of degrees of freedom per atom is exactly three. 

The intermediate phase shows up if a certain amount of medium range order (MRO) is 
allowed. This is realized in the SIC A construction by setting the quantity 64/62 non zero, 
i.e. edge-sharing (ES) tetrahedra ^^4/2 leading to four-membered rings -62^4 can now be 
created at the growing cluster steps. This means also that if stress should be created when 
f is increasing, then it should be only in ring structures and not by two corner-sharing 
connected BA4/2 tetrahedra. 

Two transitions are now obtained for every SICA step. The first one hes always around 
the mean coordination number fd — 2.4 where the number of floppy modes / vanishes. The 
second transition is located at fc2 and is a new feature of rigidity theory. When starting from 
a floppy network and progressively stiffening the network and requiring self-organization, the 
network will reach a point beyond which stressed rigid bondings outside of ring structures 
can not be avoided anymore. This point is a stress transition. We show in fig. 2 the 
I — 2 result where / = at fd = 2.4 for different fractions of ES tetrahedra, defining an 
intermediate phase Af. We should also stress that even for a non-zero ES fraction, / = is 
always obtained at f = 2.4 for / = 2. From this analysis, it appears that the first transition 
at fci does not depend on the ES fraction, as well as the fraction of stressed rigid clusters in 
the structure. In fig. 2, the probabihties of the related stressed rigid clusters for a non- zero 
ES fraction can of course be obtained from the fioppy and isostatic ones since the sum of all 
probabilities is equal to one. 

To ensure continuous deformation of the network when B atoms are added while keeping 
the sum of the probability of floppy, isostatic rigid and stressed rigid clusters equal to one, 
the probability of isostatic rigid clusters connects the isostatic solid line at fc2- Stressed 
rigid rings flrst appear in the region fd < f < fc2 while chain-like stressed clusters (whose 
probability is proportional to 63) occur only beyond the stress transition, when 63 7^ 0. This 
means that within this approach, when f is increased, stressed rigidity nucleates through the 
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network starting from rings, as ES tetrahedra or small rings. It is easy to see from fig. 2 that 
the width Af = fc2 — ^ci of the intermediate phase increases with the fraction of ES. This 
can be extended to any MRO fraction (fig. 3) and it shows that Af is almost an increasing 
function of the ES fraction as seen from the result at SICA step I — 4. Since there is only a 
small difference between allowing only four-membered rings (ES) (lower dotted line) or rings 
of all sizes (upper dotted line) in the clusters, we conclude that the ES mostly determine 
the stress transition and hence the width. 

Finally, one can see from fig. 2 and the insert of fig. 3 that the probabihty of isostatic 
clusters is maximum in the window Af, and almost equal to 1 for the even SICA steps, 
providing evidence that the molecular structure of the network in the window is almost 
stress-free. The point at which f — shifts slightly around f = 2.4 with the SICA step I 
(e.g. see the insert of fig. 3). 

B. Constraint free energy 

Prom the cluster distribution obtained by the SICA, it is also possible to compute the 
constraint- related free energy, following the approach reported by Naumis^''. Here, we have 
kept from the internal energy of the network only the part related to energy of the elastic 
deformations of the network, removing the contributions from the harmonic vibrations of 
the atoms and the anharmonic contributions which are irrelevant for our purpose (however, 
see^^). This idea is also consistent with the work of Duxbury and co-workers who showed 
that the number of fioppy modes behaves as a free energy for both rigidity and connectivity 
percolation^^ . The entropy of the network can be evaluated as a Bragg- Williams term from 
the distribution of cluster probabilities pi at step I. 

Ml 

Fi^Ui- TSi = NksTfi + NkBT J] Inp^ (3.2) 

i=l 

where fi — S — -nSp is the number of fioppy modes computed following equ. (2.7). Since 
we have expressed the latter quantity as a function of the mean coordination number f and 
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since the probabilities can also be expressed as a function of f, the free energy can be plotted 
as a function of f in the case of self-organization. Figure 4 shows the constraint related free 
energy Fi for two non-zero ES fractions. 

It appears from the figure that the second transition at =fc2 (the "stress transition" ) is 
a first order transition while the first transition at fd is weakly second order. Moreover, as 
one can see from the insert, the first transition at fci progressively becomes first order in 
character when the rate of edge-sharing tetrahedra rj is increased. On the other hand, the 
jump of Fi^^ at f = fc2 is reduced when the ES fraction is increased (insert of fig. 4). 

Both curves show a marked minimum of Fi in the range [fc2,2.667] at a certain coordi- 
nation number fg which signals an equilibrium state with respect to cross-linking. A major 
consequence of this result is that one may expect phase separation in the stressed rigid 
region leading for increasing cluster sizes to nano-scale phase separation in the network 
backbone. Close to f = 2.667, the structure should therefore be made of B-poor clusters 
having the statistics of the local Fi minimum but also compensating B-rich clusters in order 
to still satisfy equ. (2.6). This structural change is driven by the entropic term appearing in 
equ. (3.2) since the energy of the elastic deformation of the network is zero in the stressed 
rigid phase. With increasing ES fraction, this equilibrium state shifts to the value 2.667. 
For a ES fraction of 1, equilibrium state and stoichiometric composition merge together. 
Experimental evidence of nanoscale phase separation is discussed in the following. 

IV. DISCUSSION 

A. The Boolchand Intermediate phase 

As mentioned above, chalcogenide glasses are the first systems that have been carefully 
studied and the intermediate phase defined by the two transitions has been discovered by 
Boolchand in the context of self-organization^^'^^. SICA provides therefore a benchmark to 
check the results obtained. To be specific, Raman scattering^^'^^ probes elastic thresholds 
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in binary SixSei-x or GcxSei-x- The germanium or silicon corner sharing mode chain 
frequencies change with mean coordination number f of the glass network. These frequencies 
exhibit not only a change in slope at the mean coordination number fd — 2.4, but also a 
first order jump at the second transition fc2- In germanium systems, the second transition 
is located around the mean coordination number of 2.52 whereas fc2 — 2.54 in Si based 
systems. For both systems, a power-law behavior in f — fc2 is detected (see fig. 5) for f > fc2 
and the corresponding measured exponent is very close to the one obtained in numerical 
simulations of stressed rigid networks^^. Moreover, these results clearly correlate with the 
vanishing between fd and fc2 of the non-reversing heat flow AHnr (the part of the heat flow 
which is thermal history sensitive) in MDSC measurements^^'^^. 

The study of stoichiometric compounds such as SiSe2 or GeSe2 also leads to better 
understanding of medium range order. In the former, '^^Si NMR experiments have shown 
that most of the tetrahedra were part of long edge-sharing chains^" in the glass. SiSe2 has 
different crystalline polymorphs which all exhibit a strong edge-sharing tendency^^. The 
high temperature phase is made of 100% edge-sharing tetrahedral^^, while different phases 
display a distribution in terms of NMR functions (where the subscript k — 0,1,2 
refers to the number of tetrahedra shared by edges on a tetrahedron), but with a majority 
of E^"^^ structures^^'^^. In the SiSe2 glass, the fraction of E^'^^ has been found to be of 
the order^^ of 0.53. On the other hand, low temperature crystalline GeSe2 has no edge- 
sharing tetrahedraP^, but glassy GeSe^. exhibits a companion Raman line associated with 
edge-sharing tetrahedra^). 

The SICA approach has shown that the width Af of the intermediate phase increases 
mostly with the fraction of ES tetrahedra. We stress that the width should converge to 
a lower limit value of Af compared to the step 1 = 2, therefore one can observe the shift 
downwards when increasing I from 2 to 4. This limiting value is in principle attained for 
/ — > oo, or at least for much larger steps^^ than I = 4. For Si-Se, Af = 0.14 is much more 
sharply defined than for Ge-Sc (Af = 0.12) consistent with the fact that the number of ES 
tetrahedra is higher in the former 
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B. Nanoscale phase separation 



Phase separation effects in glasses exhibit usually pronounced changes in physical prop- 
erties and most studies have focused on thermally driven heterogeneities which can in some 
cases display bimodal glass transition temperatures. Here the separation effect results from 
a change in network connectivity which has its origin in the free energy minimum at fg. 

It appears that these nano-scale phase separation have been revealed from compositional 
trends^^ of the glass transition temperature Tg because they display a maximum in Tg close 
to the stoichiometric concentration. Such a feature has been observed in Ge-Se^^, As-Se^^ 
alloys, but not in Si-Se^^. These Tg trends can be compared with spectroscopic (Raman^^, 
Mossbauer^^) data which also give evidence of broken chemical order, suggesting that the 
structure of stoichiometric glasses such as GeSe2 or As2Se3 is made of a chalcogen rich 
majority phase and a compensating Ge- or As- rich phase. Furthermore, in the metal rich 
phase Ge-Ge or As-As bonds are present. 

The difference between the Si-Se and the Ge-Se glass lies in the following. Since the 
Si-Se has a much higher ES fraction compared to Ge-Se, the value of its corresponding local 
constraint free energy equilibrium will lie very close to the value f = 2.667 (see insert of 
fig. 4). The Ge counterpart will have the same minimum at a lower value in the range [fc2, 
2.667] because of fewer ES tetrahedra thus favouring the emergence of the chalcogen-rich 
phase when f is increased. Our last comments brings us back to constraint counting. Since 
Si-Se is more weakly constrained than Ge-Se due to the higher amount of ES^, the glass 
transition temperature of the stoichiometric glass will be higher compared to Ge-Se. 

V. APPLICATION TO FAST IONIC CONDUCTORS 

One interesting field of application of cluster construction and constraint counting al- 
gorithms is the field of fast ionic conductors (FIC)''^'''^, which has received considerable 
attention in the last fifteen years because of potential applications of these solid electrolytes 
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in all solid state electrochemical devices and/or miniaturized systems such as solid state 
batteries. An important step forward has been made by replacing the oxygen in usual oxide 
glasses by more polarizable chalcogenide atoms (S,Se mostly) which has increased the dc 
conductivity^^ in these systems by several orders of magnitude^^, up to a value of about 
10~^ Q-^.cm"^. Surprisingly, the extension of constraint theory from network chalcogenide 
glasses such as AsxSei-x to ionic glasses has received little attention and to our knowledge, 
has been only reported for a few oxide glasses^^'^^. Elastic percolative effects in these types 
of networks have not been studied so far with the network change in solid electrolytes, al- 
though it is certainly fundamental for the understanding of the mobile carriers' motion since 
various models of conductivity^^'^^ stress the importance of the mobility /i in the contribu- 
tion to conductivity. Obviously, the mobility is related to a local mechanical deformation of 
the network"^*^, allowing a cation to move through holes in the structure. In terms of rigidity, 
one may therefore expect that the mobility // in a stressed rigid solid electrolyte should be 
substantially lower compared to the cation mobility in a floppy one, because in the latter 
floppy modes allow a local low energy deformation. The percolative effect of mobility should 
certainly show up in this kind of network so that the conductivity a should display some 
particular behavior in the stress-free intermediate phase and at the two transitions. 

The SICA approach can be applied to the present solid electrolyte case by considering 
the simplest binary conducting glass, which is of the form (1 — x)SiX2 — XM2X with X an 
anion of Group VI (X=0,S,Se) and M an alkah cation (M=Li,Na,K,. . . ). The free carriers 
are the M® cations. The local structure in these glasses can be determined by many different 
expertiments and is usually described in terms of so-called and units, derived from 
NMR data^^. The former corresponds to the usual sihca tetrahedron made of one silicon and 
four Group VI atoms at the corner [e.g. 81864/2] while the latter has one additional anion 
bonded to the alkali cation [e.g. SiSef/2^0,^] that does not connect anymore to the network 
(flg. 6). Although it is yet not clear what is the coordination number of the alkah cation^°, 
it can be assumed that the strongest interaction of the alkali cation is the one related to the 
NBO. This means that the effective coordination number of M is taken as one, as suggested 
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by different authors^^'^^. 

Starting from the local and units (with respective probabilities 1 — p and 
p = 2x/{l — x), the SICA probabilities can be evaluated for different steps of cluster sizes 
following the procedure described previously and taking into account the 1-fold M cations^^. 
When constraint counting is performed, it appears that the creation of a — connec- 
tions leads to a stressed rigid cluster, while the Q'^ — and — connections yield 
respectively isostatically stressed and floppy clusters. The SICA results show again that 
the intermediate phase exists only when a non-zero fraction of small rings is allowed in the 
self-organized structure. The corresponding results are displayed in fig. 7 for / = 2 and 
work on higher SICA steps is in progress^^. Rigidity nucleates here in a way opposite to 
network chalcogenides. The vanishing of the number of floppy modes defines the upper limit 
of the intermediate phase, while stressed rigidity outside of ring structures disappears for 
X > Xci- This is consistent with the fact that the network is stressed rigid at low modifier 
concentration. However, the base network glass is stress free. In Si02, the Si-O-Si angle 
distribution is quite wide leading to broken bond-bending constraints on oxygen^'', while in 
the SiS2 and SiSe2 glasses, the structure is mostly made out of edge-sharing SiSe4/2 or 
SiS4/2 tetrahedra^*^ which are weakly stressed (i.e. ric < 3.667 a value which would be ex- 
pected in a dendritic network at the mean coordination number f = 2.667 or concentration 
X — 0.333) due to the 4-membered ring correction coming from the counting of redundant 
constraints^. 

Increase of the alkali content leads to an increase of fioppiness. In the oxide system 
(1 — x)Si02 — XM2O, the width Ax should be very small or zero since the fraction of ES 
tetrahedra in the oxide systems is almost zero^^. Still, percolative effects are expected at 
the concentration x = 0.2 corresponding to the transition from rigid to floppy networks, a 
transition that has been observed in sodium tellurate glasses^^. In sulfur and selenide glasses 
such as (1 — x)SiS2 — xNa2S, the width should be much broader because of the existence 
of the high amount of edge-sharing tetrahedra in the 8182 or 8i8e2 base networks^^. For 
glasses with a high amount of ES tetrahedra, the lower limit at Xd of the intermediate phase 
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is expected to decrease down to x — for the limiting case rj — 1. In the sulfur base glass, 
"^^Si NMR have shown that the fraction of ES tetrahedra should be about 0.5, slightly higher 
than in the selenide analogous system^^'^^. 

Prom fig. 7, for rj — 0.5 one should observe a window of about Ax — 0.09. Unfortunately, 
conductivity, structural and thermal results^^ on these systems are only available for an alkali 
concentration x > 0.2. However, in the different silica based glasses, a rigidity transition has 
been observed^^ at the concentration x — 0.2 which should provide guidance for forthcoming 
studies in this area. 

Finally, temperature effects should be observable close to this transition. Since the con- 
centration of alkali free carriers ul depends on the temperature (the higher the temperature, 
the higher n^), an increase of the temperature T should lead to a decrease of the number of 
network constraints, the fraction of intact bond-stretching constraints of the alkali atom 
being proportional to 1 — ul- Consequently, a shift of the mechanical threshold (/ = 0) to 
the higher concentrations should result from an increase of T. 

VI. SUMMARY AND CONCLUSIONS 

To summarize, we have shown in this article how stress change in molecular systems 
can be described using both cluster construction and constraint counting. This permits 
to go beyond the usual mean-field approach of rigidity and to obtain the two observed 
rigidity transitions. We have found that there is a single transition from floppy to rigid 
networks in a certain number of structural possibilities. An intermediate phase appears 
when a fraction of medium range order is allowed in self-organized networks and the order 
of the underlying phase transitions is first and second order, and depend also on the fraction 
of ES. Nano-scale phase separation appears in the stressed rigid phase and is driven by the 
cluster entropy. This separation leads to Group Vl-rich clusters and chalcogen rich clusters 
when the stoichiometric composition is attained. Finally, extension of this approach to ionic 
conductors has been emphasized and should motivate new developments in this field. 
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List of Figures 



Prom the short range order molecules {I — 1) yielding the mean-field result 
to all clusters at I = 2 and some MRO produced at / = 3. Isomers start to 
be created at step / = 3. Note the creation of medium range order such as 
edge-sharing tetrahedra or rings. Each boundary atom is counted half. ... 23 
Probabihty of floppy, isostatic stressed and stressed rigid clusters, as a func- 
tion of the mean coordination number f for different possibilities of medium 
range order. The solid line corresponds to the dendritic case while the broken 
lines correspond to a respective ES fraction at the stress transition of 0.156, 
0.290 and 0.818. The flUed square indicates the stress transition at the point 
fc2 and the filled circle the point fd that does not depend on the ES rate (see 
text for details). Por clarity, we have removed the probabihties of stressed 

rigid clusters for non-zero ES fractions 24 

Width of the intermediate phase as a function of the fraction of medium range 
order at the stress transition for / = 2 (sohd fine), / = 3 (dashed line) and 
I = 4 (dotted lines). The lower dotted line corresponds to a structure where 
only ES tetrahedra have been allowed. The insert shows the probabihty of 
isostatic clusters with mean coordination number for / = 3 (dashed line) and 
/ = 4 (dotted line), compared to the shaded region defined by the Af from 
SICA analysis. The point defined by / = is shifted compared to Z = 2. . . . 25 
Upper panel: Pree energy Fi of the system Ai^^Bx as a function of the mean 
coordination number f for different fractions of edge-sharing units rj. Open 
circles: rj — 0.29, filled circles: rj — 0.56. The insert shows the equilibrium 
coordination number fg with respect to the ES fraction (see text for details). 
Lower panel: the first derivative of the free energy f/^'* with respect to f as 
a function of f 26 
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Variation of the BSe4/2 (B=Ge,Si) corner-sharing mode frequency normahzed 
to 1 in Raman spectroscopy with respect to the mean coordination number 
f. The sohd vertical hnes define the intermediate phase in Ge-Se^^ while the 
lower solid hne and the dashed line define it for Si-Se, foUowing^^. The Si 

Intermediate phase is larger than the Ge one 

The two local structures and in {1 — x)SiX2 — XM2X glasses"^^, with 

X — 0,S,Se and M — Na, K, The structure has one non-bridging 

anion 

Critical concentrations in {l — x)SiX2 — xM2X glasses, with X = 0,S, Se and 
M — Na, K,. . as a function of the fraction of edge-sharing (ES) in the base 
81X2 glass. The insert shows the corresponding width of the intermediate 
phase as a function of the same quantity 
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FIG. 1. Prom the short range order molecules (Z = 1) yielding the mean-field result to all clusters 
at Z = 2 and some MRO produced at Z = 3. Isomers start to be created at step I = 3. Note the 
creation of medium range order such as edge-sharing tetrahedra or rings. Each boundary atom is 
counted half. 
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Mean coordination number r 
FIG. 2. Probability of floppy, isostatic stressed and stressed rigid clusters, as a function of 

the mean coordination number f for different possibilities of medium range order. The solid line 

corresponds to the dendritic case while the broken lines correspond to a respective ES fraction at 

the stress transition of 0.156, 0.290 and 0.818. The filled square indicates the stress transition at 

the point fc2 and the filled circle the point fci that does not depend on the ES rate (see text for 

details). For clarity, we have removed the probabilities of stressed rigid clusters for non-zero ES 

fractions. 
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Fraction of MRO at the rigidity transition 
FIG. 3. Width of the intermediate phase as a function of the fraction of medium range order at 

the stress transition for I = 2 (sohd hne), I = 3 (dashed line) and / = 4 (dotted lines). The lower 

dotted line corresponds to a structure where only ES tetrahedra have been allowed. The insert 

shows the probability of isostatic clusters with mean coordination number for Z = 3 (dashed line) 

and I = 4 (dotted line), compared to the shaded region defined by the Af from SIC A analysis. 

The point defined by / = is shifted compared to Z = 2. 
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2 2.1 2.2 2.3 2.4 2.5 2.6 

Mean coordination number r 




FIG. 4. Upper panel: Free energy Fi of the system Ai^^Bx as a function of the mean coordina- 
tion number f for different fractions of edge-sharing units i]. Open circles: rj = 0.29, filled circles: 
ij = 0.56. The insert shows the equilibrium coordination number fg with respect to the ES fraction 
(see text for details). Lower panel: the first derivative of the free energy Fi^^ with respect to f as 
a function of f. 
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2.2 2.3 2.4 2.5 2.6 
Mean coordination number r 

FIG. 5. Variation of the BSe4/2 (B=Ge,Si) corner-sharing mode frequency normahzed to 1 in 
Raman spectroscopy with respect to the mean coordination number f. The solid vertical lines 
define the intermediate phase in Ge-Se^^ while the lower solid line and the dashed line define it for 
Si-Se, foUowing^^. The Si Intermediate phase is larger than the Ge one. 
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Rate of ES in the base glass 

FIG. 7. Critical concentrations in (1 — x)SiX2 — XM2X glasses, with X = 0,S,Se and 
M = Na,K, . . ., as a function of the fraction of edge-sharing (ES) in the base 81X2 glass. The 
insert shows the corresponding width of the intermediate phase as a function of the same quantity. 



29 



